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Abstract 

A method is proposed for the calculation of diffusion constants for one-dimensional 
maps exhibiting deterministic diffusion. The procedure is based on harmonic inver- 
sion and uses a known relation between the diffusion constant and the periodic orbits 
of a map. The method is tested on an example map for which results calculated by 
different other techniques are available for comparison. 



1 Introduction 



Deterministic diffusion plays an important role in a variety of physical phe- 
nomena and applications, e.g. particle confinement in magnetic fields [1], con- 
ductivity in metals [2], etc. A general theory of diffusion for systems with 
continuous time is presented in [3]. In systems with periodic phase-space struc- 
ture, diffusion may be induced by chaotic dynamics inside the elementary cells, 
causing an irregular jumping between the cells. The investigation of determin- 
istic diffusion is an interesting task both from the theoretical point of view, 
concerning the question of how the typical stochastic properties arise from a 
purely deterministic dynamics, as well as with regard to its various applica- 
tions in different physical contexts. 

It has been observed that even the simplest model systems - one- dimensional 
maps - which exhibit deterministic diffusion are capable of explaining the be- 
haviour of relevant physical systems [4,5]. For such maps, it has been demon- 
strated that the diffusive properties are related to the periodic orbits of the 
map restricted to an elementary cell [4-7] . In particular, the diffusion constant 
characterizing the diffusive spreading can be obtained from the leading zero 
of a zeta function consisting of a product over contributions from all peri- 
odic orbits of the reduced map. Similar to the dynamical zeta functions which 
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appear in semiclassical quantization problems, this product usually does not 
converge. In former work [4-6], the convergence problem was handled by the 
application of cycle expansion techniques. However, cycle expansion strongly 
depends on the existence of a symbolic dynamics and the shadowing of long 
orbits by combinations of short ones. 

On the other hand, it has been demonstrated [8,9] that the convergence prob- 
lems of the zeta functions and corresponding response functions arising in 
semiclassical quantization problems can be circumvented by the application 
of harmonic inversion, which is a universal method in the sense that it does 
not rely on a symbolic code or other special properties of the system. The 
method has recently been extended from Hamiltonian systems to maps [10]. 
In this paper we demonstrate that harmonic inversion can also be used for the 
calculation of diffusion constants of one-dimensional maps. In analogy with 
the procedure for semiclassical quantization, we rewrite the zeta function as a 
response function, from which we construct a periodic orbit signal including a 
finite set of periodic orbits. The analysis of the signal by harmonic inversion 
yields the zeros of the underlying zeta function, from which the diffusion con- 
stant is calculated. We apply this technique to an example map for which the 
periodic orbits can easily be calculated, and which has already been examined 
with the help of cycle expansion techniques [11], rendering a comparison with 
previous results possible. 

We start by describing the general procedure. We consider one-dimensional 
maps / : IR — > R defined by an anti-symmetric function /, 

Xn+l = f{Xn) , f{-x) = -f{x) (1) 



for points x inside the elementary cell [—1/2, 1/2), and a discrete translational 
symmetry outside the elementary cell, 

f{x + n) = /(£) + n , n e Z . (2) 



Prom this map, a reduced map / : [—1/2,1/2) — > [—1/2,1/2) can be con- 
structed by restricting the dynamics to the elementary cell via the definition 



x^x -[x + l/2\ (3) 
f(x)^ fix) -[fix) + 1/2], (4) 

where [ • ] denotes the integer part of the quantity in the brackets. Each pe- 
riodic orbit of the reduced map fix), consisting of a sequence of fixed points 
{xi, X2, ■ ■ ■ , Xn}, corresponds to an orbit of the full map fix) which returns to 
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an equivalent point, 



(Jr. 



(5) 



The diffusion constant associated with the map is determined by the lattice 
translation a after n applications of the map, averaged over initial conditions 
in the elementary cell 



D = lim —{a^)n . 



(6) 



As was shown in Ref . [6] , for hyperbolic or almost hyperbolic maps, the diffu- 
sion constant can be determined from the periodic orbits of the reduced map 
(4) by considering the zeta function 



z{z,p)= nn 1- 



m=0 p 



lA I A"^ 



(7) 



where the second product runs over all primitive periodic orbits of the reduced 
map; Up and Ap = Yh=i f'i^i) topological length and the stabihty of 

the orbit, respectively, ap denotes the lattice translation of the corresponding 
orbit of the full map, and /9 is a free parameter. According to Ref. [6], the 
diffusion constant D can be expressed in terms of the leading zero Zq{(3) of the 
zeta function as a function of /?, 
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a/32 



(8) 



13=0 



Additionally, from the symmetry of the map, it follows that [dzo/dp]j3=o — 0. 
Using these relations, an expansion of In 2:0 for small /3 yields 

\nzM^-P'D, (9) 



which can directly be used to calculate the diffusion constant from the leading 
zero for small f3. 

The problem that remains is how to actually calculate the leading zero, since 
the expression (7) for the zeta function usually diverges. In past works, this 
problem has been tackled by application of cycle expansion techniques. As 
an alternative, we demonstrate in the following how the zeta function can be 
expressed in terms of a response function, and how the leading zero can then 
be obtained by harmonic inversion of a periodic orbit signal constructed from 
this response function. 
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2 Hcirmonic inversion for maps 



We start by defining the quantity w — —i \nz and rewriting the zeta function 
(7) as a function of w, 

ZK^-nn l-^TlA^ . (10) 



The density of zeros of the zeta function (10) can be written in the form 

pK/5) = --Im^(«;,/3), (11) 

TT 



where the response function g{w,(3) is given by 

g{w, (3) = go{w, Z?) + ^ In Z{w, (5) (12) 



with gQ{w,j3) a smooth, slowly varying function of w. Using the definition 

piwup pap/3 



we rewrite the second term of the response function as 



-lnZ(^,/3)=EEH^p)7fTM (14) 



g oo ^(m) 

m=0 p -1- — 

oo oo 

= EEE(-^^.)(4"^T (15) 

m=0 p r=l 
oo oo ^iwrup ^rcTp/S 

= E E E(-^M (16) 

m=0 p r=l rVI P 

piwrii,,, .,<T,,nS oo 

- E(-^^o) |. I E , (17) 

where the first sum in (17) now runs over a// periodic orbits (p.o.) including 
multiple traversals of primitive orbits, and Uq denotes the topological length 
of the underlying primitive orbit. After carrying out the sum over m, the 
response function finally reads 

g{w, (5) = go{w, (5) + E(-^^o) _ . , . (18) 

po I "-po I 
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On the other hand, the response function can be written in terms of the zeros 
Wk of the zeta function (10) and their multiphcities d^, 



^K/3) = E (19) 



Like the zeta function (10), the expression (18) for the response function will 
usually not converge in those regions of the complex variable w, where the 
poles of g{w, j3) are located, and the zeros of the zeta function cannot directly 
be calculated from (18). Instead, we apply harmonic inversion to obtain an 
analytical continuation [8,9]. The poles Wk{P) of the analytical continuation of 
g{w, P) in (18) are the zeros of the analytical continuation of the zeta function 
Z{w,P) in (10). To apply harmonic inversion, we now construct a periodic 
orbit signal by Fourier transformation of the response function (18), 

C{s,P) = — J g{w,P)e-'^^dw. (20) 



Considering only the second, oscillating part of g{w,(3), this results in 

C{s,P)^Y: I S{s-n,^). (21) 



po I ^ '-po I 



The corresponding signal resulting from Fourier transformation of Eq. (19) 
reads 

C{s,P)^-i^dke-^^^^^^\ (22) 



which contains the zeros Wk{l3) of the zeta function (10) as frequencies. The 
zeros Wk{(3) can now be determined by adjusting the periodic orbit signal (21) 
- including all orbits up to a maximum topological length rimax - to the form 
of the corresponding signal (22) by harmonic inversion. According to Eq. (9) 
and with the relation w = —i \nz, the diffusion constant can be obtained 
from the leading frequency wq using the relation 

D^^ wo{(3) , (23) 



which for small (3 should be independent of (3. 
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3 Results 



As an example we consider the following map, which has also been examined 
in Ref. [11]: For points x inside the interval [—1/2, 1/2] the map is defined by 

/(£) = £(1 + 2|2£|") , (24) 

where a is a free parameter with the restriction a > — 1. Outside the elemen- 
tary cell, the map is continued according to Eq. (2). Using Eqs. (3) and (4), a 
corresponding reduced map /(x) can be constructed. 

As was discussed in Rcf. [11], the periodic orbits of the reduced map f{x) can 
be described by a complete symbolic dynamics with a three-letter alphabet. 
For — 1 < a < the dynamics is completely hyperbolic. The orbit (corre- 
sponding to the fixed point x = 0) is infinitely unstable (i.e., Ap = oo) and 
does not contribute to the response function (18). For a > 0, the orbit be- 
comes marginally stable with Ap = 1, and the dynamics becomes intermittent. 
Furthermore, it was shown in [11] that at a = 1 the system undergoes a phase 
transition from normal to anomalous diffusion, and the diffusion constant as 
defined by Eq. (6) equals for a > 1. 

For a > 0, the marginally stable orbit would give a diverging contribution 
to the zeta function (7) as well as to the periodic orbit signal (21). As was 
discussed in Ref. [7] , it is appropriate to simply omit this orbit from the zeta 
function, as was also done in Ref. [11]. In Ref. [11], the diffusion constant of 
the map (24) as a function of a was determined using explicit cycle expansion 
expressions for the average in Eq. (6) rather than making use of Eq. (8). In 
the following, we will determine the diffusion constant by harmonic inversion 
as outlined above. 

We calculated the diffusion constant for different values of the parameter a in 
the region —0.5 < a <2 using the periodic orbits up to topological length = 
13 (but omitting the orbit). For each value of a, we performed calculations 
for different small values of (3 and checked the convergence of the results for 
the diffusion constant for /3 — > 0. The results for the diffusion constant are 
presented in Figure 1. In the region a < 0.5, our results are in excellent 
agreement with those from Ref. [11] obtained with the periodic orbits up to 
topological length A^ = 10. For a > 0.5, the different techniques applied 
in Ref. [11] yielded quantitatively different values for the diffusion constant; 
here, our results show a good qualitative agreement with the different cycle 
expansion results. This is even true when we used the reduced set of orbits 
up to topological length A^ = 10 in our calculations: The function D{a) in 
Fig. 1 becomes slightly more flat, but still shows a good qualitative agreement 
with the results of Ref. [11]. As in [11], it was not possible to reproduce the 
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theoretical value L> = for o; > 1, but one only obtains an asymptotical 
convergence to this value. 

In our calculations, we observed that the condition that for given a the ex- 
pression (23) for the diffusion constant should become independent of /3 for 
/5 ^ was well fulfilled for a < 0. However, with increasing a > 0, the results 
became numerically unstable if (3 was too close to zero, so that we had to 
extrapolate from the results for (3 around 0.5. The convergence of the frequen- 
cies and amplitudes in the harmonic inversion procedure itself was also very 
good for a < 0, but convergence was increasingly hard to obtain for a > 0. A 
criterion for the reliability of the results for a given value of a is the closeness 
of the leading frequency for /5 = to its theoretical value = 0, as well as 
the closeness of the corresponding amplitude do (see Eq. (22)) obtained by 
harmonic inversion to its theoretical value 1 (the multiplicity of the leading 
zero of the zeta function). The absolute value of the leading frequency for 
/3 = and the corresponding amplitude as a function of a are presented in 
Figure 2. The absolute value of the leading frequency was smallest for « < 
but stayed below \wq\ ~ 10~^ for all values of a. The corresponding ampli- 
tudes were nearly exactly equal to their theoretical value 1 for a < 0, but 
began to deviate from their theoretical value for a > 0.3, indicating that the 
harmonic inversion results became less well converged. A possible reason for 
the deviations may be that the density of frequencies becomes too large for 
the signal length chosen. It can be expected that the accuracy of the results 
will be improved by including longer orbits. 

In Rcf. [11] it was proposed to truncate the periodic orbit sum by introducing 
a stability cutoff and to use the stability as an ordering parameter for the 
cycle expansion. It is worth noting that the harmonic inversion procedure has 
the advantage that no discrete ordering parameter is needed at all. However, a 
stability cutoff can also be used with the harmonic inversion technique because 
very unstable orbits give only a small contribution to the signal which can be 
neglected. 



4 Conclusion 

In conclusion, we have demonstrated how the diffusion constant of one-dimen- 
sional diffusive maps can be determined from a finite set of periodic orbits us- 
ing harmonic inversion techniques. Starting from a zeta function whose leading 
zero possesses a known relation to the diffusion constant, we have constructed 
a periodic orbit signal which was analysed by harmonic inversion. The dif- 
fusion constant could then be determined from the leading frequency of the 
signal. Since, in contrast to other methods such as cycle expansion, harmonic 
inversion does not depend on any special properties of the map, the general 
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procedure should work in the same way for aU chaotic maps. Wc have tested 
the method for a simple map whose periodic orbits can easily be determined 
and for which results for the diffusion constant obtained by cycle expansion 
techniques are available in the literature. Our results are in excellent agree- 
ment with the ones from cycle expansion. 
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Fig. 1. Results for the diffusion constant of the map (24) as a function of the 
parameter a, calculated by harmonic inversion of the periodic orbit signal (21) 
including all periodic orbits up to topological length 13 (with the exception of the 
orbit). 
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Fig. 2. Test of the convergence of the harmonic inversion procedure: Absolute 
value of the leading frequency Wq for /3 = and the corresponding amplitudes do 
obtained by harmonic inversion, as a function of the parameter a. The closeness 
of these values to their theoretical values ■wo{P = 0) = and do{f3 = 0) = 1 are a 
measure for the reliability of the results for the diffusion constant obtained for the 
respective value of a. 
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